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ABSTRACT 

We analyze the thermal conductivity of ions (equivalent to the conductivity of phonons 
in crystalline matter) in a neutron star envelope. 

We calculate the ion/phonon thermal conductivity in a crystal of atomic nuclei 
using variational formalism and performing momentum-space integration by Monte 
Carlo method. We take into account phonon-phonon and phonon-electron scatter- 
ing mechanisms and show that phonon-electron scattering dominates at not too low 
densities. We extract the ion thermal conductivity in ion liquid or gas from literature. 

Numerical values of the ion/phonon conductivity are approximated by analytical 
expressions, valid for T > 10^ K and 10^ g cm~^ ^ P ^ 10^^ g cm^"^. Typical magnetic 
fields B ^ 10^^ G in neutron star envelopes do not affect this conductivity although 
they strongly reduce the electron thermal conductivity across the magnetic field. The 
ion thermal conductivity remains much smaller than the electron conductivity along 
the magnetic field. However, in the outer neutron star envelope it can be larger than 
the electron conductivity across the field, that is important for heat transport across 
magnetic field lines in cooling neutron stars. The ion conductivity can greatly reduce 
the anisotropy of heat conduction in outer envelopes of magnetized neutron stars. 

Key words: dense matter - stars: neutron 



1 INTRODUCTION 



Neutron stars with strong magnetic fields B > 10^^ G 
are expected to have an anisotropic surface tempera- 
ture distribution (see, e.g., Gepert. Kueker fc Page! |2004| . 



l2006l : iPerez-Azon'n. MiraUes J.A. fc Pons J.A.I bOOef i. It is 

thought to result from the anisotropy of thermal conduction 
in a neutron star envelope and manifests in periodic modula- 
tion of X-r ay emission observed from some spinning ne utron 
stars (e.g.. lBurwitz et al.|[20o3 : iHaberlllioOTl : IHoI|2007| ). The 
thermal conductivity k± across the magnetic field can be 
muc h smaller than the conductivity k h along the field (see, 
e.g.. lYakovlev fc Kaminkerlll994l : |Potekhin 1999. 1. Both con- 
ductivities ar e crucial for rnodeling of cooling magnetized 
neutron stars (iGep pert et al."2006l:|P erez-Azorm et al ]|2006l . 
IPage. Geppert fcW cbcr 2006, and references therein). 

In a non-magnetized envelope, th ermal energy is mainly - 
transported by elec trons (k ~ k^; iFlowers fc Itohl 1 19761 : 
iPotekhin et al.lll999l ). The thermal conductivity of ions, Ki, 
in which we include the ion conductivity of an ion gas or liq- 
uid and the phonon conductivity of a crystalline ion solid, 
is usually much smaller than Kc- However, a strong enough 
magnetic field B can easily suppress Kc±, making Ki the 
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leading thermal conducti vity across the magnetic field (e.g., 
IPerez-Azorm et akllioO^ ). 

In this paper we analyze different regimes of the ion 
thermal conductivity in a neutron star envelope. Physical 
conditions in the envelope are described in Sect. [J] The ther- 
mal conductivity in the magnetized envelope is outlined in 
Sect. |3l where we give a sketch of anisotropic electron heat 
conduction and nearly isotropic ion/phonon heat transport. 
The ion conductivity due to ion-ion (phonon-phonon) scat- 
tering is studied in Sect, [l] The ion conductivity owing to 
ion-electron (phonon-electron) scattering is analyzed in Sect. 
[H In Sect. |6] we estimate the values of magnetic fields which 
can strongly affect ion heat transport. Section [7] contains a 
discussion of our results. Conclusion is presented in Sect. 13 



2 COULOMB PLASMA IN A NEUTRON STAR 
ENVELOPE 

We consider a model of a neutron star envelope 
whic h combines several models of ground state matter 
( Ncge le fc Vautherinlll973: lOvamatsul 1993i and others) and 



which is summarized bv lHaensel. Potekhin fc Yakovle vll2007l 
(their Appendix B). Typical densities and temperatures of 
practical interest are indicated in Fig. [l] and outlined below 
in this section. The shaded regions I-V are domains, where 
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Table 1 . Ion conduction regimes I-V in Fig. [T] 



Regime 


State of ions 


Leading 






scattering* 


I 


Weakly coupled ions 


ii 


II 


Classical ion liquid 


ii 


III 


Low-T liquid, high-T crystal 


ii (phph) 


IV 


Quantum ion crystal 


ie (phe) 


V 


Very cold quantum crystal 


ie (phe) and 






possibly others 



'Symbols i, e, and ph refer to ions, electrons, and phonons, 
respectively. 



the ion thermal conductivity has different character. These 
domains are summarized in Table [T] and described further 
in Sects. |3H5] The last column of Table [T] gives the leading 
ion (phonon) scattering mechanism. In addition, in Fig.[T]we 
show the boundaries of the regions, where ion heat transport 
is strongly affected by the magnetic fields Bm = 10^'^ G and 
lO^** G (the regions below and to the left of the lines) . 

At any density p the plasma is assumed to consist of 
electrons and single species of fully ionized atoms (bare 
atomic nuclei). In the inner neutron star envelope, at den- 
sities higher the neutron drip density pnd ~ 4.3 x 10^^ 
g cm~^ [shown by the v ertical dot-dashed line in Fig. [1] 
iNegele fc Vautherin|[l973l ] the matter contains also free neu- 
trons. 

W e describe the p aram eters of dense matter following 
§2.1 in iHaensel et al.l (|2007l ). The charge neutrality of the 
matter implies 



rii = Tic/Z 



(1) 



where Z is the charge number of the nuclei, while Wo and rii 
are the electron and ion number densities, respectively. The 
mass density of the matter can be estimated as 



; riiA'mu 



(2) 



where ruu is the atomic mass unit and A' is the number of 
nucleons per one nucleus. For p < pnd, one has A' = A, 
where A is the number of nucleons confined in one nucleus. 
For p > Pnd, one has A' = A + A" , where A" is the number 
of free (unbound) neutrons per one nucleus. The ion mass is 



mi = A niu ■ 



(3) 



First, we describe the properties of electrons (Sect. tnj) 
and ions (Sect. 12. 2p neglecting the effects of the magnetic 
fields, and then (Sect. [Sj we outline the magnetic effects. 



2.1 Electrons 

Under the conditions of study (Fig.[T]), the electrons are typi- 
cally strongly degenerate. The electron Fermi momentum pp 
and Fermi wavenumber fcp are 

PF = fifcp = ft {S-K^ric)^^^ . (4) 
The electron relativity parameter can be written as 

1/3 



= ^ ^ 1.009 ( 



(5) 



o 




6 8 10 12 

logio P (g cm-3) 

Figure 1. (color online) Temperature-density diagram for the 
ground-state neutro n star crust (smooth-composition model, 
iHaensel et"al] [looj); Ti is the temperature for the onset of 
Coulomb coupling of ions into ion liquid; Tm is the crystalliza- 
tion temperature of ions; Tp is the ion plasma temperature; the 
vertical line shows the neutron drip density. Shaded are regions 
I-V of different regimes (Table[T]l of the ion thermal conductivity. 
Also, we show the boundaries of the regions where ion heat trans- 
port can be strongly affected by the magnetic fields = 10^"^ G 
and lO'^* G. See text for details. 



where pe = p/10^ g cm ^ . The electron Fermi energy and 
the effective mass at the Fermi surface take the form 



ep = mcC (1 + x^) 



2U/2 



^ mc (1 + a:: 



2\l/2 



ep 
^2 



and the electron Fermi velocity reads 



The electron degeneracy temperature is 



Tp = 



ep 



fc: 



"^"^-=5.93xl0« 



1-^x2-1 K 



(6) 
(7) 

(8) 



where fce is the Boltzmann constant. The electrons are 
strongly degenerate as long as T ^ Tp. 

Finally, the electron Thomas- Fermi screening wavenum- 
ber (inverse plasma screening length of Coulomb interaction 
due to polarizability of strongly degenerate electrons) is 



fcTP 



1/2 



fcp , 



(9) 



where a is the fine structure constant. 



2.2 Ions 

The properties of the classical ion plasma are determined by 
the Coulomb coupling parameter 

(10) 

where 
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Ion thermal conductivity 3 



(47rni/3)-'/^ 



(11) 



is the ion sphere radius. The plasma ions form a nearly ideal 
gas as long as F ^ 1 which corresponds to T ^ Ti = 
Z'^ / {aiks) ■ At lower T they transform smoothly into a 
strongly coupled ion liquid. 

It is also important to define the ion plasma frequency 

LOp={'i-nniZ'^e^/miY^'^ , (12) 
and the associated ion plasma temperature 



hiOp 



7.832 X 10" 



£1 ^ 
A' A 



1/2 



(13) 



displayed in Fig. [T] Quantum effects in ion motion become 
very important for T <C Tp. 

The conditions of i on crystallization are described, e.g., 
in iHaensel et al.l |20o3), §2.3.4. A classical Coulomb liquid 
of ions crystallizes at the temperature 

where Fm « 175. Large zero-point vibrations of ions at very 
high densities can reduce Tm or even prevent crystallization. 
Actually, this effect can be important only for a plasma of 
light nuclei (hydrogen and helium) which we do not consider 
here. Ions are usually assumed to form body-centered cubic 
lattice, but, in fac t, the face c entered cubic lattice is also 
possible (see, e.g.. lBaikoll2002l ). We study phonon conduc- 
tivities for both lattice types and find, that they are nearly 
identical. The same is true for other kinetic and thermo- 
dynamic p roperties of such crystals (see, e.g., iBaiko et al.l 
Il998l . l200ll ). 



3 THERMAL CONDUCTIVITY IN A 
MAGNETIZED PLASMA 

As outlined in Sect. [T] we are interested in the thermal 
conductivity of magnetized neutron star envelopes. A mag- 
netic field B makes thermal conduction anisotropic (e.g., 
[Yakovlcv & Kaminkcr 1994; Potokhin 1999). The total ther- 
mal conductivity can be written as 

K = Kc -|- Ki, (15) 

where Kc and Ki are the conductivities of electrons and ions, 
respectively. In the inner envelope of a neutron star one 
should also add the thermal conductivity of neutrons. How- 
ever, it is largely unexplored and will be neglected here. 
Generally, both conductivities, Kc and Ki, are anisotropic 
tensor quantities. Nevertheless, a typical neutron star mag- 
netic field B ~ 10^^ - 10" G c an strongly affect the electron 
thermal condu ction (see, e.g., lYakovlev fc Kaminkeilll994l : 
|Potekhin|[l999l ) but weakly affects the ion thermal conduc- 
tion (see Sec. Therefore, we will take into account the 
effect of the magnetic field on Kc but neglect its effect on ki. 
Let us emphasize that our aim is to study the ion conduc- 
tivity. We describe the electron conductivity here only for 
comparison with the ion one. 



3.1 Thermal conductivity of electrons 

An anisotropic electron conduction in a magnetic field B 
is characterized by the thermal conductivity k^ii along the 



magnetic field, by the conductivity Ke± across the field, 
and by the Hall conductivity KeH (perpendicular to B and 
to the temperature gradient). The effects of the magnetic 
fields on electron thermal conduction are twofold (see, e.g., 
lYakovlev fc Kaminkej|l994 IVentura fc Potekhin|[200l] ') . 

First, there are classical effects associated with elec- 
tron Larmor rotation about magnetic field lines. Their ef- 
ficiency is characterize d by the electron ma gnetization pa- 
rameter ojgTe (see e.g. lYakovlev fc Kaminker 19941. where 
Uf, — eB/(m*c) is the electron gyrofrequency, and Te is the 
electron relaxation time at B = 0. In this classical approxi- 
mation 



— ^1 ^ — 



3m* 



AO) 



Ko± = 7—7^ , (16) 

and Kcu = Ka±ijJgTe. If the electrons are strongly magne- 
tized {ujgTc ^ 1), their fast Larmor rotation greatly reduces 
Ke± and KcH. In this limit, KcH becomes a non-dissipative 
quantity (independent of Tc). 

Second, electron transport can be modified by quan- 
tum effects associated with the structure of electron Lan- 
dau levels in the B-field. These effects are especially pro- 
nounced at sufficiently low densities p < pbc, at which 
degenerate electrons populate the only one (ground-state) 
Landau level. The critical density psc can be estimated 
as pBe = 7.045 X 10^ {A/Z){B g cm'^ (see, e.g., 
|Potekhin| [l999). At higher densities p > pBe, the electrons 
populate other Landau levels and, when the density in- 
creases, all conductivity coefficients, especially Ke|| and KeX, 
oscillate with growing p around their classi cal values (1161) i n 
response to the population of new levels (|Potekhinl 119991 ). 
These quantum oscillations are typically not too strong and 
can be neglected here for our semi-quantitative considera- 
tion of electron thermal transport. For simplicity, we will 
restrict ourselves to high densities p > pBc and neglect the 
quantum effects, but retain much stronger classical effects. 

In order to evaluate the electron conductivity co- 
efficients in our approximation we need the expression 
for Te. We will calculate i t including electron -ion scat- 
tering (iPotekhin et al ] |l999l : iGnedin et al.|[200]] ) [although 
neglecting freezing of U mklapp processes d iscussed by 
iRaikh fc YakovlevI Il982l : iGnedin et all l200lt and a re- 
cent ly revised contribution from electron-electron scatter- 
ing (|Shternin fc YakovlevI [200^ ). For this purpose, we will 
use the Mat thiessen rule 
Chapter 7 in lZimanlll960t ). where 
five relaxation times (at B = 0). Strictly speaking, such an 
inclusion of electron-electron collisions into transport coeffi- 
cients (|16|) in a magnetized plasma is approximate. Similarly, 
at T < Tp electron-ion collisions become essentially inelas- 
tic and Eq. (|16p in a magnetized plasma is quantitatively 
inaccurate even if electron-electron collisions are neglected. 
However, we employ these approximations for illustrating 
the importance of the ion thermal conductivity. 



T"ei* +r^e^ (see, e.g., §10 of 
Tci and Toe are partial effec- 



3.2 Thermal conductivity of ions 

As already mentioned above, we will neglect the effects of 
magnetic fields on the thermal conductivity of ions, and de- 
scribe ion transport by one coefficient of thermal conduc- 
tivity, Ki (for 5=0). Estimates of magnetic field strengths 
which make ion heat transport anisotropic are given in Sect. 
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[6l Their typical values are ~ 10^'' G, much higher than the 
magnetic fields which induce a strong anisotropy of the elec- 
tron transport. 

The ion thermal conductivity can be presented as 



«:i = (^ii + '^ie''") 



(17) 



where kh and partial ion conductivities due to ion-ion 

and ion-electron collisions, respectively. In fact, Kio appears 
to be important only for crystallized ions (see Sec. [7]). There 
are several regimes of ion conduction realized at different T 
and p (domains I-V in Fig. [T] summarized in Table [TJ . We 
analyze them in subsequent sections. 

Let us remark, that should not be confused with 
Kei, the electron thermal conductivity due to electron- 
ion (electron-phon on) scattering; Kci is well known (e.g. 
iGnedin et allbOOll ) and gives main contribution to the elec- 
tron transport. 



4 ION CONDUCTION DUE TO ION-ION 
COLLISIONS 

4.1 Gaseous phase 

If r < 1 (r > Ti in Fig. [l]), the ions form a nearly ideal 
Boltzmann gas and the ion thermal conductivity can be cal- 
culated from the standard Boltzmann equation taking into 
account ion-ion and ion-electron Coulomb collisions. In this 
case, ion-electro n collisions are negligible and, according to 
iBraginskil (|l963D . 



«: i = 3.9 



4ko r-'/'/Aii 



(18) 



where kq = oJ-p rii O; is a convenient normalization con- 
stant, 



3 m,'/' (fee Tf'^ 



(19) 



" 47ri/2 (^e)4 Aj. 

is the effective relaxation time due to ion-ion collisions, = 
In [l/(V3r'''')] is the appropriate Coulomb logarithm, and 
Oi is the ion-sphere radius given by Eq. (|lip . In this regime, 
ion conduction is realized through weak Coulomb collisions 
of almost free ions. In the limit of F <C 1 the Coulomb 
logarithm is sufficiently large reflecting long-range nature of 
Coulomb interactions in the weakly coupled plasma. For a 
moderate coupling (F ~ 1), one has An ~ 1 because of the 
onset of strong ion screening. 



4.2 Classical Coulomb liquid 

Here we outline Kii in a classical strongly coupled ion liq- 
uid (1 J; F < Fm), where the ions are no longer free 
but are mostly confined in Coulomb potential wells. These 
results are expected to be especially suitable in domain 
II in Fig. [T] In this case kJ? can be calculated using the 
molecular dynamics formalism. The physics of ion trans- 
port becomes essentially different from the gaseous case (see, 
e.g., Bernu fc Vieillefosselll978l : iPierleonT Ciccotti fc BermJ 
1 19871 '). A nearly fi-ee ion motion in the gas is replaced by oscil- 
lations in Coulomb potential wells, with occasional hopping 



from one well to another (see, e.g., lDaligaultll2006l ). Hop- 
ping transitions can contribute to the effective thermal con- 
ductivity along with thermal energy transfers via Coulomb 
interaction of ions vibrating in neighboring potential wells. 
Such heat transport processes can be described in terms of 
ion-ion scattering via emission and absorption of phonons 
in the ion liquid; they c an be studied by molecular dy nam- 
ics technique (see, e.g.. iMcGaug hev fc Kavianvll200a ). The 
resulting thermal conductivity can be written as 



0.4 Ko, 



(20) 



where k, is a dimensionless, sl owly variable function of F . 
Thi s function was calcula ted by [Bernu fc VieiUefossl (|l978D 
and lPierleoni et al.1 1| 19871 ) for F ~1, 10, and 100, and the re- 
sults are in a good agreement (see Fig. [2]). For our simplified 
semiquantitative analysis, it is sufficient to set this function 
equal to a typical value k« = 0.4. It gives us a tempera- 
ture independent thermal conductivity nlj which combines 
smoothly with the conductivity p8|l in the gaseous phase 
at F ~ 1. Expressing in the form of familiar estimate 
Kii ~ niTiifc|r/mi [see Eq. (|18|) ]. we obtain an estimate 
of the effective relaxation time, t^^ ~ r/tjp, which gives 
at F ~ 1. 



4.3 Coulomb crystal 

Now let us study kh in domains III-V in Fig.[T] which mainly 
refer to crystalline matter. In this case, ion transport can be 
described using the formalism of phonons (ele mentary ex- 
citations of the crystalline lattice; IZimaniri960l . Chapter 1) 
and Kii is determined by phonon-phonon scattering (absorp- 
tions and emissions of phonons, ph, which can be schemat- 
ically presented as ph 4- ph ^ ph and ph ph + ph) . For 
the parameters of dense matter we are interested in (Fig.[T]), 
the approx imation of almost uniform electron b ackground 
holds (e.g. JPotekhin. Chabrier fc Yakovlev|[l997l '). Then the 
phonon spectrum and anharmonic terms (which determine 
phonon-phonon scattering) of Coulomb crystals can be com- 
puted with high precision (e.g.. lCohen fc Keffeill 19551 : iDubinI 
1990). Accordingly, one can accurately calculate kh but this 
calculation is complicated and we restrict ourselves by semi- 
quantitative estimates. Note that phonon scattering in the 
crystal is similar to that in the ion liquid although the prop- 
erties of phonons in the li quid and the c rystal can be some- 
w hat different (e.g., jj71 in Landau fc Lifshitz 1993i . and §24 
in lLifshitz fc Pitaevskilll98ol )! 

For a classical Coulomb crysta l (T^^/ S ^ T < Tm) we 
employ a simple estimate bv lZimaJ l|l96Cf ) [his Eq. (8.2.16) 
in Chapter 8] 



ci latt rnim Cs 



72 r 



27 T 



(21) 



Here, 7 — 9 In Tp/(? In — 0.5 is the Griineisen constant and 
we assume that the typical phase velocity of phonons is Cs ~ 
ciJp/(3gBz) [we replace the B rillouin zone by a sphere w ith 
the radius qsz ~ (67r' ni)''^^]. iPerez-Azorfn et al. I l|2006t ) es- 
timated Kii^'^^^ from the same expression but using 7 = 2, 
characteristic for terrestrial solids (where the electron back- 
ground is essentially non-uniform). Introducing our conduc- 
tivity normalization constant kq, from Eq. (|2H) we obtain 
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Ion thermal conductivity 



2437r 



1/3 



(22) 



The thermal conductivity Hu in a crystal is solely deter- 
mined by Umklapp phonon-phonon scattering (the Peierls 
theorem; e.g., Chapter 8 in Ziman 1960), in which the sum 
of two (either initial or final) phonon wavevectors goes out- 
side the first Brillouin zone. In addition, there are so called 
normal scattering processes, where the sum of wavevectors 
remains in the first Brillouin zone. With decreasing T (or 
increasing p) the Umklapp processes become frozen out 
because they involve phonons with uj ^ uo^ (whereas the 
amount of such phonons at T ^ Tp - in the quantum, low- 
temperature limit - is exponentially small). The freezing 
strongly enhances kh. 

Calculations of Kii at T <C Tp are difficult; some es- 
timates of K jj unde r diffe rent assumptions are presented in 
Chapter 8 of lZimanI (119601 '). However, for a neutron star crust 
it is sufficient to use a simpler estimate, consisting in mul- 
tiplying the conductivity (|22p for a classical crystal by an 
appropriate enhancement factor, 



exp(/?rp/r). 



(23) 



where /3Tp ~ Tp/3 corresponds to the estimation of the 
lowest frequency of phonons participating in Umklapp pro- 
cesses. We will use this expression everywhere in domains 
III-V. The enhancement factor exp (/JTp/T) can be inaccu- 
rate at T ^ Tp but the main contribution into the total 
ion conductivity ki at these low T comes from ion-electron 
scattering (Sect. [5]), so that the inaccuracy does not strongly 
affect the total ion conductivity ki. 

Note that the phonon therm al conductivity can be pre- 
sented in the form (|Zimanlll960l . §1 in Chapter 7) 



(24) 



where Ci is th e phonon (dimens ionless) heat capacity per one 
ion (see, e.e:.. lBaiko et al.|[200lf ): Cg ■p^ ii'p/(3gBz) is a typical 
phonon group velocity, and Lph is an effective phonon mean 
free path. Equation (|23|l gives 



^ph 



: 20 ai 



exp(/3Tp/r). 



(25) 



We remark that d « 3 for T > TJ?, an d d ~ {T/Tpf 
for T < Tp/3 (see, e.g.. iBaiko et alj|200ll ). Let us empha- 
size that the effective mean free path Lph should be calcu- 
lated neglecting normal phonon-phonon scatte ring processes 
which do not affect phonon heat transport jZimanl Il960l . 
Chapter 8). 

At T < Tp, an important contribution into can re- 
sult from phonon-impurity scattering. At very low tempera- 
tures, phonon scattering by crystal boundaries can become 
significant (the neutron star crust can be a polycrystal made 
of small monocrystals). However, in view of our ignorance 
concerning the actual nature and distribution of impurities 
and polycrystal structures in the crust, these mechanisms 
will not be considered in the present paper. In the high- 
temperature limit {T > T'p/3), we obtain a not too large 
mean free path Lph ~ 20ai. It indicates that the presence of 
not very abundant impurities (a few percent by number) or 
bulky monocrystals will not affect the ion thermal conduc- 
tivity. In the low-temperature quantum limit (T <C Tp/3) 
the phonon mean free path is restricted by electron scatter- 
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Figure 2. (color online) Dimensionless ion-ion thermal conduc- 
tivity Kii/fto for a classical ion plasma (T 3> Tp) as a func- 
tion of r. Circles sho w results of equilibrium mo lecular dynamics 
(EMD) simulations of lBernu &: Vieillefossell978l (BV 78), squares 
and diamonds present equilibrium and non equilibrium molecu- 
lar dynamics results of IPierleoni et al.l lligSTl l (PCB87 EMD and 
NEMD), respectively. The solid line is our interpolation, Eq. I I26I I. 
Dotted, dot-dashed and dashed lines are the conductivities in do- 
mains I-III, given by Eqs. JTSj, I I20I I and I I22I I. respectively. 



ing (Sec.[5l and can be estimated as L ~ 200ai [Sec. 15.31 Eq. 
(|43[l ]. In this case heat is mainly conducted by long wave- 
length phon ons (Sec. | 5.2I> . whose scattering by impurities 
is inefficient (jZimanlll96'o[ Chapter 8, §3). This means that 
the impurities are not very important at low temperatures 
as well. 



4.4 Interpolation expression for Kii 

To facilitate the implication of the above results we can sug- 
gest the following interpolation formula for kh in all domains 
I-V in Fig. m 



Kii — 



+ K 



1/2 



Aco [v-'/\n' [2 + l/(x/3r='''')] 



+ + [v/nf exp (2/3rp/r)} 



1/2 



(26) 



with K* « 0.4 and P ~ 1/3. Here, we assume a smooth 
temperature variation of Kii, without any jump or break at 
the melting point T = T^. Our assumption is made in anal- 
ogy with the electron thermal c onductivity which does not 
show any peculiarity at T = Tm (|Baiko et al.lll99^ ) because 
of the importance of multi-phonon electron-phonon scatter- 
ing processes in crystalline lattice at T < Tm and because 
of incipient ion- ion correlations in ion liquid at T > Tm. 
Another argument in favor of a relatively smooth behavior 
of the conductivity near the melting point is given by the 
presence of shear modes in a strongly coupled Coulomb liq- 
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uid (|Schmidt et al.lll997h . that is typical for a crystal rather 
than for a liquid. 

Figure [2] shows the dependence of a dimensionless ther- 
mal conductivity kh/kq in a classical ion plasma on the 
Coulomb coupling parame ter F. Circles are the results of 
iBernu fc Vieillefosse (1 19781) (BV78) . Squares and diamonds 
are calculations of Pierleoni et al.l (|l987l ) carried out us- 
ing equilibrium molecular dynamics and nonequilibrium 
molecular dynamics techniques (PCB87 EMD and PCB87 
NEMD), respectively. Dotted, dot-dashed and dashed lines 
present the thermal conductivities calculated, respectively, 
from Eqs. (|18p . (|20|l and (|22p . which are valid in domains I, 
II and III. Our interpolation ((26}, plot ted by the solid line, 
is in g ood a greement with t h e resu lts of lBernu fc Vieillefosse! 
l|l978t ) and iPierleoni et al] (119871 ). In particular, the inter- 
polation predicts a mini mum of kh/kq at F ^ 7, w hich 
agrees with the estimate of lBernu fc Vieillefosse! (|l978l ) that 
the minimum takes place in the vicinity of F = 10. The 
largest, but still acceptab le deviation of the interpo lated 
ku/kq from th e resu lts of iBernu fc Vieillefosse! l|l978! ) and 
IPierleoni eTsU (Il987l ) occurs at F ~ 100. It can be explained 
by a significant contribution of multi-phonon processes [ne- 
glected in Eq. (|22|l ]. as well as by an inaccuracy of Eq. (I22|) . 
which is just an order-of-magnitude estimate of the ther- 
mal conductivity of ions in the crystalline phase that we use 
to determine the thermal conductivity at F = 100 in our 
interpolation. 



5 PHONON CONDUCTION DUE TO 
PHONON-ELECTRON SCATTERING 

Now let us discuss kic, which is the partial ion thermal con- 
ductivity due to ion-electron s cattering in Eq. (|17ll. Thi s 
process was not considered by IPerez-Azon'n et al.! (|2006! ). 
As already mentioned above, this mechanism appears to be 
important only for the crystalline lattice at T < Tm, and 
becomes negligible at higher T as compared to ion-ion scat- 
tering. 



5.1 General formalism 

We will mainly focus on the case in which T < Tm and 
the phonon formalism is appropriate, so that can be 
treated as the phonon thermal conductivity produced due 
to phonon-electron scattering (due to absorptions and emis- 
sions of phonons by electrons, ph-|-e^e and e— >ph-|-e). We 
expect that the formalism described below is valid for cal- 
culating Kie in temperature-density domains III-V displayed 
in Fig. [T] and listed in Table [T] The conductivity be 
calculated using the va riationa l method, described, for in- 
stance, in Chapter 7 of IZimanI (|l960! ). The variational ex- 
pression reads 



5c, 



(27) 



In this case, v = (Q, s) enumerates phonon modes, Q 
is a phonon wavevector, and s is a polarization index; 
a phonon frequency will be denoted by Uv. Furthermore, 
= l/[exp(fia;„/feBr) — 1] is the phonon equilibrium occu- 
pation number, and v,^ = diOv/dQ is the phonon group ve- 



locity; S'coii is the entropy generation rate in phonon-electron 
scattering. A variational function $i/ describes a weak devi- 
ation of the phonon distribution function from the equilib- 
rium distribution n^. The variational estimate of Kio given 
by Eq. (|27|) reaches maximum for the exact solution of the 
phonon transport equation. For phonon-electron sca ttering, 
the entropy generation rate can be expressed as (|Ziman! 
1 1960! . Chapter 8, §9) 



'S'ci 



j ^QpQ.fcdQdfcdfc'. 



(28) 



For the sake of simplicity, polarization indices are omitted 
here. Wavevectors k and k' refer to electrons. In calculating 
Kie the electron distributions are assumed to be equilibrium 
ones. Finally, Vq^^ is the differential transition probability, 

calculated in equilibrium. 

The derivation of Kie is simplified l|Zimanlll960! . Chapter 
8, §9) by the equality of the squared matrix elements for the 
phonon absorption and emission processes. Employing the 
simplest suitable variational function = Q ■ u [u being 
the unit vector along the temperature gradient), one gets 



■iZe 



where 



pIvf 



i-jT Z^e^ riiAe 



(29) 



(30) 



are auxiliary quantities, which are formally similar (but not 
identical) to the electron electric conductivity and the elec- 
tron effective relaxation time due to electron-photon scat- 
tering, and 



Ac = 



2p| 



1- ^ 



1 
2kF 



An 

exp[-M/(g)] 



exp(- 



[1 - exp(-'2:^)]^ 



(31) 



is formally similar to the Coulomb logarithm for the elec- 
tron transport [although it differs from the real Coulomb 
logarithm that contai ns instead of Q^; see, e.g 



iBaiko fc Yakovlev!ll995! |. Since a and are not real elec- 
tric conductivity and relaxation time, kic and a do not 
obey the Wiedemann-Franz law. The integration in (I31|l is 
done over positions of both electrons on the Fermi surface, 
with wavevectors fc and fc' (dSlfe and dSlfe/ are respective 
solid angle elements). Furthermore, q = fc — fc' is an elec- 
tron momentum transfer in a scattering event, and Q is 
the phonon wavevector that is equal to q reduce d to the 
first Brillouin zone (see Chapter 1 of IZimanlll960! for more 
details). Specifically, Q = q — g, where g is a reciprocal 
lattice vector which realizes such a reduction (absorbs a 
momentum excess). Let us remind that the processes with 
g 7^ belong to Umklapp phonon-electron scattering pro- 
cesses, while the processes with g = are normal phonon- 
electron processes. The quantity |0((?)| ~ F{q)/{q^ ~\- /ctp) 
in Eq. (|31|l is a Fourier transform of the screened Coulomb 
electron-ion interaction; /ctf is the Thomas-Fermi screening 
wavenumber given by Eq. ((9)1; F{q) is the nucleus form fac- 
tor which accounts for the proton charge distribution within 
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the atomi c nucleus; Wjq) is the do ubled Debye- Waller factor 
(see, e.g.. lBaiko fc YakovlevllT995h . In addition, we have de- 
fined z,j — huJs(Q)/k-BT and the phonon polarization vector 
e-u = es{Q). 

We have performed extensive computations of Ac and 
/tie from Eqs. (|29p - pi|l by the Monte Carlo technique. In 
these calculations, the positions of electrons (fe and fe') on 
the Fermi surface are randomly selected at every step in long 
Monte Carlo runs. For every selection, we calculate q, a re- 
spective inverse lattice vector g, and a phonon wavevector 
Q. The phonon eigenfrequencies uja{Q) and polarization vec- 
tors es{Q) have been determined then for a Coulomb crystal 
of ions immersed in the rigid electron background (see, e.g., 
iBaiko et al.ll200lh . Small (but finite) widths of energy gaps 
in electron energy spectrum at the intersections of the elec- 
tron Fermi surface with boundaries of Brillouin zones have 
been neglected. We have calculated for body-centered 
cubic and face-centered cubic crystals, and the results turn 
ou t to be the same. T he same conclusions have been reached 
bv lBaiko et al.l (|l998l ) with regard to the electron transport 
coefficients determined by electron-phonon scattering. 

When the temperature decreases, the number of 
phonons, which can efficiently participate in phonon- 
electron scattering, becomes smaller. They can only be the 
phonons with small frequencies [huJa{Q) ~ ksT] and, hence, 
with small wavevectors Q near the center of the Brillouin 
zone. The regions on the Fermi surface, which mainly con- 
tribute to the integral H3ip . become narrower, and the accu- 
racy of our direct Monte Carlo calculations gets lower. At 
sufficiently low T direct Monte Carlo calculation is inefficient 
and we replace it with a semi-analytic consideration (Sect. 
15. 2[) . In Fig. [T] the T — p domains, where direct Monte Carlo 
calculations are effective, are labeled as III and IV, while 
domain V will be described semi-analytically. We combine 
Monte Carlo and semi-analytic results and produce a useful 
practical fit expression in Sect. 15.31 



5.2 Umklapp processes at low temperatures 

Let us analyze phonon-electron Umklapp processes {g 7^ 0) 
at low temperatures (domain V in Fig. [IJ . In this case it is 
sufficient to consider small phonon wavenumbers Q <C gez, 
and we can rewrite Eq. (|31|l in the form 



A.^ = 



pI 



8 7r2 rtiikjiT 



E 



F\9) 



exp[ 



9<2feF 



2 



2fcF/ 



2 exp(-0^) 



[1 — exp(— 2^ 



(32) 
T2' (33) 



where the sum is over all non-zero reciprocal lattice vectors 
inside the Fermi sphere {g < 2fcF). Using symmetry prop- 
erties, one can show that [g ■ e^)^ in the integrand can be 
replaced by 3^/3. 

At low temperatures we are interested in, the main con- 
tribution into Ig comes from those positions of k and k' 
on the Fermi sphere, for which k — k' ~ g. This case can 
be stu died semianalytically following IPethick fc ThorssonI 
j 19971 ) who considered analogous problem in the neutrino 
emission due to electron-nucleus collisions. Let us introduce 
a coordinate system whose z axis is directed along g. Let 



us further define k by its polar and azimuthal angles 
and if, and fc' by its polar and azimuthal angles and 
ip' , and introduce /j, = cosi9, /i' = cosi?'. The allowed po- 
sitions of fc should concentrate on a ring {0 ^ (p < 27r) 
with ~ /io = g/(2fcF); while the respective positions of fc' 
should concentrate on a complementary ring (tp' f^ip = ifo, 
i9' ~ TT — i9). Let fco be a position of fc exactly on the ring, 
and fc') = fco — g be a respective {(p' — <p) position of fc' on 
the complementary ring. Small variations of fc and fc' around 
these positions can be defined by = fJ,~fJ,o, A^' = fi'—fig, 
Aip' — Lp' — ifo, Ak — fc — fco, and Afc' = k' — kg. Using 
cylindrical coordinates, we get then 



Afc 



Akj 



Afc'|«fcFA/i', Afcl ; 

AkL 



Mo 



MO 



fcp Afi, 



1 — /Iq kpAifi' . 



(34) 



Because the phonon wavevector is Q = Afc — Afc' , we obtain 



Q\\ 
Q± 



fcp {Afi' + Afi) , 



MO 



Mo 



fcp (A/i' - A/i) 



1 — Mo A(f' . 



(35) 



The Jacobian of the transformation from variables 
(A/i, A/i', A(p') to (Q||, Q±, Qy,) is 1/2 Mofcp- It is then easy 
to pass from coordinates {Q\\, Q±, Qip) to spherical coordi- 
nates of the vector Q. Because only small-Q phonons are 
available in the cold crystal, we can safely extend the inte- 
gration over Q to infinity. As a result, we obtain 

f2n 



F JO 

exp( 



d^c 



dgg" 



(36) 



[1 - exp(-z^)]^ 

At Q <C gaz in the Coulomb crystal we have two nearly 
acoustic phonon modes (21.2 ~ fici^2{Q) Q/ksT, Q be- 
ing a unit vector along Q) and one nearly optical mode 
(23 ~ ftajp/fcaT 3> 1) that does not contribute to Ig. The 
integration over dQ for acoustic modes gives 



f 

Jo 



exp(-2i,2) 
[1 - exp(-2i,2)]' 



■dQ = 4!C(4) 



ksT 

fiCi,2(Q) 



15 



fcsT 



_ftCi,2(<?)_ 

and the integral over gives 27r. Therefore, 



(37) 



4 = 



where 



87V^ 

45 

647r^ 
45 



-T 



ksT 

h 

kBT 



dS^c 



-1/5 



(38) 



(39) 



is an average phonon group speed that can be written as 
= Auo-p/qw- The calculation yields A ~ 0.36 for both 
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bcc and fee Coulomb crystals. Finally, we have 



45 mi Ljp \ kp 



1 

g 

2kp 



9<2fcF 

exphW^C?)]. (40) 



As in the ease of electron transport (see, e.g. JPotekhin et alj 
1 19991 ). normal phonon-electron processes give small contri- 
bution into scattering rate and can be neglected. The same 
is true for higher temperatures. 

Notice in passing, that in the same manner we can ob- 
tain an asymptotic form of the Coulomb logarithm for the 
electron electrical conductivity due to electron-phonon scat- 
tering, 

-2 „2 _ -, / ^ N 2 



A? 



7!" ^tez 9bz 1 I T 



9 mio^p ftp A'' 



2kF 



3<2kp 



1- ^ 



9 
2kF 



eM~W{q)]. 



(41) 



However, in this case an effective phonon group velocity is 
Ca' = A' uJp/qBz with A' ~ 0.40 for both bcc and fee lattices. 
The difference in the values of A and A' results from aver- 
aging different powers of Cs(Q). Our asympto tic expression 
for A n agrees with previous calculations (e.g. jGnedin et aU 
120011 ). 

5.3 Interpolation expression 

In order to approximate the thermal conductivity kio in do- 
mains III-V we present it in the familiar form 

Kic = CimCsI/ph, (42) 

with Cs = (^p / {3qBz) ■ We have extracted the mean free path 
Lph from our numerical calculations and obtained that it can 
be described by a slowly variable function ofO — Tp/T which 
we denote by F{6). Its density dependence has been approx- 
imated by replacing the sum over reciprocal lattice vectors g 
in the asymptotic expression (|40p with the integral over the 
vectors g uniformly distributed within the spherical layer of 
the inner radius gaz and the outer radius 2fcF. Then, the 
Coulomb logarithm A^ has been adjusted to the mean free 
path Lpij using Eqs. ((29)) . (|30j, and (|42)) . While integrating, 
we have used the Debye- Waller factor which corresponds to 
a given temperature. In analogy with the approxim a tion o f 
the electron thermal conductivity by ICnedin et al.l l|200lh . 
the effect of the nucleus form factor has been taken into 
account by introducing into the integrand an additional fac- 
tor exp (uiform (?^/4 fcp) . As a result, the mean free path has 
been approximated as 

3 mi cijp p| lip F{9) 



Li 



ph 



where 



2 Z^e'im* (jBZ nc Aphe 

320 ai 26 F{e) f Apa 
(1 -Ha;?) Aphe 'Z OOT \ ~W 



1/2 



2 A, 



ph c 



F{e) 



0.014 



E^{w) 
0.03 



exp(6l/5) + 1' 



(43) 

(44) 
(45) 



w — wuw + Inform, and y 
quantity 



gBz/(2fcF) = (4^)-!/^ The 



[0.5?t-i exp(-9.1/6l) -f W-2/6'] (46) 

mea sures the efficien cy of the Debye- Waller factor 
(Ba iko fc Yakovlevl ll995l: Wform determines t he importance 
of th e form factor of atomic nuclei (see, e.g., ICnedin et al.l 
I2OOII ): u_i ~ 2.8 and u- 2 ~ 13.0 are dimensionle ss phonon 
frequency moments (e.g.. iPollock &: Hansen|[l973l ): Ei{x) = 
exp{—y) dy is the standard exponential integral 
(who se approximations are given by Abramowitz & Stcgu3 

Now let us employ the ground-state composition of 
erustal matter and the corresponding smooth composi- 
tion model to describe proton charge distribution within 
at omic nuclei in the crust (see Sect. [2] and Appendix B 
of iHaensel et al.ll2007^ ■ The radial dependence of the pro- 
ton number density within the proton core of the nucleus 
(r < _Rp) is np{r) oc [l — (r/_Rp)*p] , where _Rp and tp are 
density dependent parameters. In order to describe the effect 
of the nuclear form factor it is sufficient to set 



WCovm = 43 Xi, 



(47) 



where 



2 

•^nuc 



1 - 15/(5 + tp) + 15/(5 + 2tp) - 5/(5 + 3tp) 
a? 1 - 9/(3 -htp)+ 9/(3 -I- 2tp)- 1/(1 +tp) 



(48) 



iHaensel et al.|[2007l . Appendix B). Notice that the nucleus 
form factor and the associated parameter lOform are im- 
portant only in inner crust, where the phonon conductiv- 
ity is rather insignificant. In the outer crust, one can set 

Wioim = 0. 

We have calculated on a dense grid of T and p values 
(logjg T [K] from 5 to 9, with the step of 0.2; logj^Q p [g em~'^] 
from 5 to 14, with the step of 0.1) restricted by the condi- 
tion r 30. The root mean square relative fit error is 8%, 
and the maximum error « 25% takes place at T = 6 x 10^ 
K and p = 4 x lO^'^ gem"''. Had we restricted ourselves to 
the outer crust (p ^ 4.3 x 10 ^ gcm"'^), where the phonon 
thermal conductivity can be really significant, then the root 
mean square relative error would be 7%, and the maxi- 
mum error of 21% would occur at T = 4 x 10^ K and 
p = 2.5 X 10^^ gcm"^. For V < ~ 175 (the shaded 
region in Fig. [HJ ions form a strongly non-ideal Coulomb 
liquid, and the crystal-phonon description is, strictly speak- 
ing, invalid. Nevertheless, as far as transport properties are 
concerned, a strongly coupled Coulomb liqui d has much in 
common with a Coulomb crystal (see, e.g. , ISchmidt et al.l 
ll997l . lBaiko et al.lll998l . iHaensel et al.ll2007l ). Therefore, we 
have extended calculations (and included them in the fit) 
to r = 30, because we expect that for 30 < F < 175 the 
adopted formalism still gives correct order-of-magnitude es- 
timates of the phonon-electron thermal conductivity in the 
ion liquid (which is the thermal conductivity of strongly 
coupled ions). Let us remind, that the electron transport 
coefficients beh ave smoothly near the melting point (e.g., 
IPotekhin et al.lil999. ). 

In fact, near the melting point (at F < 175) phonon- 
electron scattering becomes unimportant for ion conduction 
(see Fig. [1] and Table [l| and an uncertainty of our approxi- 
mation of Kic does not affect the accuracy of the evaluation 
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of the total ion conductivity Ki. A formal extension of our 
fit expression (|43|) to low F < 30 leads to an unrealistic 
exponential freezing of the ion-electron scattering (an expo- 
nential growth of Kic), which does not affect Ki determined 
by ion-ion scattering. 

Let us remark that if we neglect the Debye- Waller factor 
(wdw = 0) and the nuclear form factor (uiform = 0), then 
we obtain 

(^)-4('-lt)- <«' 

Accordingly, Aph c can be treated as a familiar Coulomb log- 
arithm (which takes into account a suppression of backscat- 
tering for relativistic electrons) with the minimum electron 
momentum transfer of qbz and the maximum momentum 
transfer of 2 ftp . 



6 ANISOTROPY OF ION HEAT 
CONDUCTION 

Let us estimate magnetic field strengths which make ion 
heat transport anisotropic. In a weakly coupled Coulomb 
plasma (F <^ 1) the effect of the magneti c field on the io n 
thermal conductivity is well known (e.g.. [Braginski|[l963l ). 
Specifically, the magnetic field does not affect noticeably 
the ion transport as long as the ion magnetization parame- 
ter is small, wsiTi < 1, where uiBi = ZeB / (niic) is the ion 
cyclotron frequency. 

In a Coulomb crystal or liquid, we expect that the mag- 
netic field strongly influences ion heat transport provided it 
distorts the spectrum of heat carrying phonons. In a clas- 
sical crystal (T > Tp/S) the heat is carried by all phonons 
whose spectrum is distorted when ujbi > oJp (|Baikdl2000l : 
iHaensel et~al]|2007l . § 4.1.6b). The latter inequality wiU be 
treated as a condition for an anisotropic ion thermal con- 
ductivity; we will use it for ion liquid and solid. In the 
quantum crystal {T <^ Tp) heat is mainly transported by 
low-frequency phonons (Sec. [5^ . with u ~ lot ~ kBT/h. 
Their spectrum is affected by the magnetic field at lubi ^ i-ot 
l|Baikdl2OO0l : lHaensel et al.ll2007l . §4.1.6b). Accordingly, the 
criterion for anisotropic ion conduction can be written as 
wsi ^ uJpT/Tp. The same criteria can be applied for describ- 
ing the effects of the magnetic field on the thermodynamic 
properties of Coulomb crystal. They are in a qual itative 
agreement with the resul ts of accurate calculations (jBaikd 
l2000l : lHaensel et alllioOTl . §4.1. 6b, their Fig. 4.3), which can 
serve as an additional argument in favor of our criterion. 

Basing on the above estimates we can introduce a char- 
acteristic magnetic field Bm, which starts to affect ion heat 
transport and makes this transport anisotropic. At B < Bm 
one can neglect the effect of the magnetic field on ki and use 
the results of preceding sections. We can suggest a simple 
estimate of B^; it is valid in a wide range of densities and 
temperatures including all values of p and T shown in Fig. 

m 

X [Stt F^VAh + 1 + {Tp/Tf] '^''^ G, (50) 
where rii is given by Eq. (|19|) with the Coulomb loga- 



rithm Aii = In 1^2 + l/(\/3F^''2)J . The magnetic field B > 

10" pI^'^ G, for which (j->Bi ^ i^p, strongly affects ion heat 
transport at any temperature. At lower field B, the condi- 
tion B ~ Bni is realized at sufficiently low or sufficiently high 
temperatures. In the case of low temperatures, a magnetized 
Coulomb crystal possesses very soft phonon modes which 
are responsible fo r heat transport and most sen sitive to 
the magnetic field (|Baikdl200Gl : iHaensel et al.ll2007l . §4. 1.6b, 
Fig. 4.3). In the high temperature plasma, the ions become 
weakly coupled and their relaxation time (I19|l rapidly in- 
creases with the growth of T, increasing a typical rotation 
angle of ions between successive collisions, and hence the 
magnetic field effect. 

The boundaries of domains, where B > Bm and Ki is 
anisotropic, are shown in Fig. [1] for Bm ~ 10^^ and Bm ~ 
10^'* G. The domains themselves are situated below and to 
the left of these boundaries. A not very strong field B — 
3 X 10^^ G (employed in Figs. |3] and [5ll does not affect ion 
conduction at all T and p in Fig. [T] 

Let us notice that at rather low densities p < psc 
(see Sec. I3.1|l the plasma electrons occupy only one or sev- 
eral Landau levels. This circumstance can strongly modify 
electron-phonon scattering. We have neglected this effect 
in the present publication. However, at these low densities 
the ion thermal conductivity (at not too low temperatures 
T > lO^K) is mainly determined by ion-ion scattering (see 
Fig.[T} and electron-phonon scattering is relatively unimpor- 
tant. At high densities the electrons populate many Landau 
levels so that the Landau level structure can be neglected 
as confirmed by calculations of electron transport properties 
(|Potekhinlll999l ). 



7 DISCUSSION 

Let us compare the ion thermal conductivity (calculated ne- 
glecting the effects of the magnetic field, as discussed in Sect. 
[3| with the electron thermal conductivities (|16p along and 
across the magnetic field B. 

Figure [3] shows the density dependence of the partial 
ion conductivities kh and Kie, the total ion conductivity 
Ki, and also of the electron conductivities along and across 
B = 3x 10^^ G, and Ko±, at three values of T = lO'^, 10^ 
and 10* K. One sees that the longitudinal electron thermal 
conductivity always dominates over the thermal conduc- 
tivity of ions. Therefore, the ion conductivity cannot con- 
tribute into heat transport along magnetic field lines under 
typical conditions in neutron star envelopes. Nevertheless, it 
can compete with Kc± and be the dominant thermal conduc- 
tivity across the magnetic field lines, especially at not too 
high densities (in the outer neutron star envelope) and tem- 
peratures. For instance, at B = 3 x 10^^ G and T = 10** K 
(the upper panel in Fig.|3)) ki does not dominate in the trans- 
verse conduction at all, but at T = lO'^"^ K (the middle and 
bottom panels) it dominates in the outer layer at densities 
p < 10® g cm~^. This density range is most important in 
the neutron star physics because it belongs to the neutron- 
star heat blanketing envelope which shields warm neutron 
star interior from the efficient cooling via thermal conduc- 
tion to the surface and t hen through the the rmal emission 
from the surface (see, e.g. . IHaensel et al]l2007l ). Higher mag- 
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Figure 3. (Color online) Thermal conductivity of the ground- 
state neutron star crust for T = 10®, lO'', and 10* K. The solid 
line is the total ion (phonon) conductivity k\. The dash and dot 
lines show k^^ and respectively. The dot-dash lines are the 
electron thermal conductivities, across and along the magnetic 
field B = 3 X 10^2 G, calculated using a simplified model of 
Sect. 13.11 The vertical dotted line shows the neutron drip density 
Pnd = 4.3 X 10^^ g cm~^. The shaded zone in the upper panel 
corresponds to a strongly coupled Coulomb liquid (T > Tm, see 
Fig. [ij. At T = 10^ and 10"^ K the ions form a crystal for all 
displayed densities. 



netic fields would stronger suppress Kai_ and widen the den- 
sity range where ion thermal conduction across B dominates 
over electron one. 

It is important to emphasiz e the efficiency of phonon - 
electron scattering (neglected by IPerez-Azorm et al.|[2"006h . 
This scattering is most efficient and determines ki at suf- 



!''''! 




I ■ ... I 
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Figure 4. (color online) Characteristic magnetic field Bi, which 
makes the ion thermal conductivity k; larger than the transverse 
electron thermal conductivity Kg^, and characteristic field Sm, 
which makes ion thermal conduction anisotropic, versus density, 
for five values of the temperature with logjg T\K\ = 6.0, 6.5, 7.0, 
7.5 and 8.0. 



ficiently high densities and low temperatures (in domains 
IV and V in Fig. [TJ, while phonon- phonon (ion-ion) scat- 
tering overtakes ion transport at lower densities and higher 
temperatures (in domains I~III) . Phonon-electron scattering 
reduces the total ion conductivity k\ in comparison with the 
conductivity that is solely determined by phonon-phonon 
scattering. In particular, for B = 3 x lO" G and T = 10* K 
(the upper panel in Fig. |3J this suppression is very efficient 
/5 > 3x 10^" g cm~'^. When the star cools, phonon-electron 
scattering becomes more important, and the layer of its dom- 
inance expands to lower p. For instance, at T = lO'^ K (the 
bottom panel of Fig. [3)l it dominates at all densities (see also 

Fig. DO). 

If we fix p and T and increase B then the heat transport 
across B will be determined by ions (ki > Kox) after B ex- 
ceeds some value B\. Figure |4] shows the density dependence 
of -Bi. At low p, this density dependence is weak but _Bi is 
a strong function of the temperature (noticeably decreases 
with decreasing T). At large p, the field _Bi becomes almost 
independent of T (all curves in Fig. U merge into nearly one 
curve) but increases with growing p. For instance, one needs 



G for the ion 
10^ g cm"^ (for 



an ordinary pulsar magnetic field B > 10 
conduction to dominate over at p 
T < 10^ K). However, the neutron star must have a typical 
magnetar field B > 10^^ G to reach similar dominance at 
p ~ 10^^ g cm^^ (for all temperatures in Fig. Also, in 
Fig. |4]we show the values of the magnetic fields, Bm, which 
make ion heat transport anisotropic for the same tempera- 
tures T. At not too high densities, p < 10^° g cm"'', the ions 
become the leading heat carriers across the magnetic fields 
B < Bm- However, at p ~ 10^^ g cm~^ the magnetic field 
B > 10^* G, which makes the transverse ion heat transport 
dominant, becomes larger than Bm if the temperatures falls 
below T ~ 10® K. However, even in that case it would be 
preferable to use the ion thermal conductivity at B = than 
to neglect the ion thermal conductivity at all. 

A magnetic field introduces a considerable anisotropy 
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loglop [g cm" ] 

Figure 5. (color online) Anisotropy of the thermal conductivity 
Kx/'^ll versus density for the magnetic field B = 3 X 10^^ G and 
three values of the temperature logj^Q T[K]=6.0, 7.0, and 8.0. See 
text for details. 

in thermal conduction of the neutron star envelope. For 
illustration, in Fig. [5] we plot the density dependence of 
the anisotropy, which is the ratio of the transverse con- 
ductivity K± = Hi + to the longitudinal conductivity 
K|| — Ki + « K^w of the matter. We take the same mag- 
netic field B = 3 X 10^^ G and the same three values of the 
temperature log^Q T—6, 7, and 8 K, as in Fig. [S] The lines 
marked 'e-|-i' are plotted taking into account both, electron 
and ion, conductivities. We see that the anisotropy can be 
very strong, and it increases with lowering T, when the star 
cools. For instance, we have k±/k^\ ~ 0.1 for T = 10* K 
and p ~ 10* g cm"^ but ~ 10"® for T = lO'' K and 

the same p ~ 10* g cm~^. To emphasize the importance of 
the ion conduction, the fines marked as 'only e' in Fig. [5] 
show the anisotropy Ke±/tie\\ produced by electron conduc- 
tion alone. At T = 10* K the ion conductivity is small and 
does not affect heat transport as discussed above (curves 
'e-|-i' and 'only e' merge). However, at lower T ion conduc- 
tion becomes important. It is clearly seen that Ki greatly 
reduces the heat conduction anisotropy at sufficiently low 
p, where ion conduction dominates in the transverse heat 
transport (see Fig.|3]and the discussion above). For instance, 
taking T = lO" K and p = 10*^ g cm and neglecting the 
ion contribution, we would have a very strong heat transport 
anisotropy k,c±/k.^^\ ~ 10~®, while the inclusion of the ion 
conductivity gives kj_ /fi:|l ~ 10"^ This means that the heat 
transport becomes much more isotropic (by approximately 
three orders of magnitude). This effect may strongly affect 
the temperature distribution in the neutron star heat blan- 
keting envelope and over the surface, the relation between 
the internal neutron star temperature and surface thermal 
luminosity, and the neutron star cooling. 



8 CONCLUSIONS 

We have calculated the thermal conductivity of ions ki 
in large ranges of the density and temperature in a neu- 
tron star envelope. We have taken into account ion-ion and 
ion-electron scattering (equivalent to phonon-phonon and 



phonon-electron scattering at sufficiently low temperatures) 
and analyzed various ion-conduction regimes. We have per- 
formed extensive Monte Carlo calculations of kic due to 
phonon-electron scattering. All our results are approximated 
by simple analytic expressions. In this way we have got a 
systematic description of all ion conduction regimes in neu- 
tron star envelopes. We have compared calculated values of 
Ki with the values of the electron thermal conductivity Kc 
which is usually the only thermal conductivity taken into 
account in neutron star envelopes (except for neutron star 
atmospheres, where heat transport is radiative). 
Our main conclusions are: 

(i) The ion thermal conductivity is much less affected 
by neutron star magnetic fields than the electron thermal 
conductivity (and we have neglected the effects of magnetic 
fields on Ki). 

(ii) The ion thermal conductivity is typically much lower 
than the electron thermal conductivity along magnetic field 
lines in the neutron star envelope. Heat conduction along 
the magnetic field lines is mostly provided by electrons. 

(iii) The ion thermal conductivity can be higher than the 
electron thermal conductivity across magnetic field lines in 
the outer envelope of a neutron star at not too high tem- 
peratures and densities (T < 10* K and p < 10^ g cm~^ for 
B ~ 3 X 10^2 G). 

(iv) The conductivity Ki at low temperatures and high 
densities is determined by electron- ion collisions (e.g., at p > 
3x 10* g cm"^ for T ~ 10^ K) although at lower p and higher 
T it is determined by ion-ion collisions. 

(v) The inclusion of Ki can strongly (by a few orders of 
magnitude) reduce large anisotropy of thermal conduction 
in magnetized neutron star envelopes, affecting thus temper- 
ature distribution in the heat blanketing envelope and over 
the stellar surface, and cooling of magnetized neutron stars. 

Our results ext end those obtained by 
IPerez-Azorfn et all (120061 '). First, we have included the 
contribution of electron-phonon scattering. Second, we have 
improved the consideration of phonon-phonon scattering. 
Finally, we have analyzed different ion conduction regimes. 
Nevertheless, we stress that our results are semi-quantitative 
and can be elaborated further. In particular, it would be 
interesting to consider phonon scattering in impure and 
imperfect crystals (by studying phonon-impurity scattering 
neglected here) , to take into account electron band structure 
effects in phonon-electron scattering, to describe accurately 
phonon-phonon scattering [without using the estimate Eq. 
(|21|) ]. and the effects of magnetic fields on the ion thermal 
conductivity. All these problems open a new and interesting 
field of neutron star kinetics which goes far beyond the 
scope of the present paper. We expect to deal with them in 
future publications. 
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